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Abstract —We demonstrate sequential mass inference of a 
suspended bag of milk powder from simulated measurements of 
the vertical force component at the pivot while the bag is being 
filled. We compare the predictions of various sequential Inference 
methods both with and without a physics model to capture the 
system dynamics. We find that non-augmented and augmented- 
state unscented Kalman filters (UKFs) in conjunction with a 
physics model of a pendulum of varying mass and length provide 
rapid and accurate predictions of the milk powder mass as a 
function of time. The UKFs outperform the other method tested 
- a particle filter. Moreover, inference methods which incorporate 
a physics model outperform equivalent algorithms which do not. 


1. Introduction 


Sequential inference has been used for measurement and 
control in many contexts including vehicle navigation ITI . 
target tracking and chemical process plant control 12. 
Sequential inference algorithms take a time-series of noisy 
measurements of a system, and produce increasingly accurate 
estimates of the system parameters or state variables in the 
form of a posterior distribution. We adapt sequential inference 
algorithms for use in automated weighing systems, where the 
system under consideration exhibits dynamics according to 
physical laws. Such methods provide dynamically updated 
estimates of mass, along with estimates of its uncertainty, 
which will be useful for the control of automated weighing 
systems. 

We use the example of a milk powder bagging system, 
where the bag is suspended from a point at its top from which 
it may swing under the influence of gravity. The bag is gradu¬ 
ally filled with milk powder while measurements of the vertical 
force component are made at discrete times. We test sequential 
inference algorithms by simulating such measurements, taking 
into account both process and measurement noise, and run 
the algorithms on the simulated data-set. We analyse the 
predictions for accuracy, precision and speed of inference, 
and draw conclusions about the most suitable algorithm. We 
stress the benefits of considering the underlying physics when 
designing sensors and control systems. 
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Fig. 1. (a) Schematic of the system. The box is a model of the bag of 

length L, which pivots at the indicated point at its top. The cross-sectional 
area is A (not indicated on the figure). The milk powder is represented by the 
shaded region, and has density p, and its centre of mass (COM) is indicated. 
The system behaves as a pendulum of length I (from pivot to COM), (b) 
Simulated data for the observed vertical component of force at the pivot during 
bag filling. 


II. Methods 


A. Physics model 

We model the bag-filling system as a pendulum with mass 
increasing at rate m (as the bag Alls), and with effective length 
I - L - Xcom, where Xcom = tnli'2-pA), p is the density of the 
milk powder and A is the cross-sectional area of the bag. The 
rate m is assumed to vary randomly (see below). Figure [T] (a) 
shows a schematic of the system. The equations of motion for 
the pendulum are simply: 

6^ to, (1) 

d>--gsm0ll, (2) 

where 0 is the angle of the pendulum, and the effective 
pendulum length / is not constant but decreases as the milk 
powder’s centre of mass rises. The system is observed via the 
vertical component of the force on the measuring device from 
which the bag is attached, which is given by: 

F — m cos 0 {itx? + g cos e) . (3) 

We simulate the system at discrete times to generate some 
example data. We propagate the dynamical variables between 
each time according to Eqs. Q and Q, and at each timestep 
evolve the flow of mass into the bag according to: 

log(OT„/lkg s^') = log(OT„_i/lkg s"') -H d„, (4) 
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where ~ N(0, 2,„) is a process noise term, which models the 
tendency of the filling-rate to slow down and speed up, whilst 
remaining strictly positive. Then each force measurement is 
simulated as; 

Ffi — Vrt, (5) 

where the measurement noise v„ ~ N 


B. Inference methods 

To compare different algorithms’ performance in estimating 
the mass of milk powder in a bag during filling we use 
three different sequential inference algorithms, each either in 


conjunction with the physics model described in Sec. II-A 


or without such a model. Each inference model is fed the 
same set of measurements {y„) (in this case equal to the 
simulated forces Ff), and is given an initial state estimate 
in the form of a prior distribution. When using the physics 
model, the algorithms estimate a state vector containing the 
components 6, oj, log(m/lkg), log(m/lkg s“'), log(L/lm) 
and log (pA/1kg m“*). Otherwise, the algorithms estimate 
log (m/1kg) and log (m/1kg s“') only. Note that log variables 
of some quantities are estimated to ensure strictly positive 
estimates of those quantities. 

Each algorithm assumes the state vector propagates between 
timesteps using the forward map f : x„ i-> x„+i. When the 
physics model is used, this map integrates the equations of 
motion [Eqs. Q and (j^], and evolves the mass as dictated by 
hi. Otherwise, only the mass is evolved. Each algorithm also 
maps a state to a measurement estimate by an observation map 
g ; X i-> y, which is given by Eq. ([^ when the physics model 
is used, and F - mg otherwise. 

We briefly describe the functioning of each algorithm below, 
and provide detailed description in Eigs. 2]|4 


1) Kalman filters: Eirstly, we use two sequential-inference 
algorithms based on the Kalman Eilter, namely the augmented- 
state Unscented Kalman Eilter (UKE), developed in Ref. IT), 
a, and the non-augmented UKE, described in Ref. a. These 
algorithms are extensions of the Kalman Eilter designed to 
perform well with nonlinear forward maps and/or nonlin¬ 
ear observation maps, on the principle that it is easier to 
approximate the probability distribution than to approximate 
(linearise) the nonlinear functions. The ability to perform well 
using nonlinear maps is vital for the current problem, since the 
equations of motion of the physics model are nonlinear. Even if 
the physics model is not used, the use of log variables to ensure 
positive estimates of the mass and mass-flow rate requires the 
forward map for these variables also to be nonlinear. 

Both algorithms approximate the state distribution by spe¬ 
cially chosen ‘sigma points’, which capture at least the first 
two moments of the distribution. The sigma points are run 
through the forward map and the transformed mean and 
covariance are used in order to perform a usual Kalman Filter 
update H). 

The augmented-state UKE differs from the non-augmented 
UKE in 3 ways: in the augmented-state UKE the sigma-point 
states are augmented with parameters representing the process 


Initialise Set the prior mean //q and covariance Ko with a 
well-motivated estimate, along with parameters a, fi and 
y which determine the distribution of sigma-points 14) . 
Provide the estimated process and measurement 
covariances: Y^d and Zy. 


for n - \...nt - 1 do 

1) Calculate sigma points and weights 




+ (l - a' +4 


N^ + A _ 

x^'^ = Bn-x + sj{Nx + 4) K„_i, 

- V(A^x + 4 )K„_i, 
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2(Nx + Ay 

for i - l...Nx, where A - {Nx + k), and Nx is the 
state dimension. 

2) Transform sigma points using the forward map: 


s-(0 


-f(x®)- 


3) Calculate the predicted mean and covariance: 






i=0 


K„ = 2, + £ (x® - (x® - iu„f. 


1=0 


4) Recalculate the sigma points using the predicted 
mean and covariance: 

X®’ = fn 


x*'* ^ fin + 




(Nx + 4) K„, 

^ fin - yf(- 


j(Nx + 4) K«. 

5) Apply the observation model to each new sigma 
point : 

y® = g(x,^iO- 

6) Calculate the predicted observation; 

INr 




end 


1=0 

7) Calculate the innovation covariance: 

2Nx 

Sn^i-y + Y (y® - y") (y® - y")^- 

/=0 

8) Calculate the cross covariance: 

2N^ 

Kf = ^W®(x®-/i„)(x®-/l„f. 

1=0 

9) Perform a usual Kalman Filter update : 
Bn — Bn "t" 

K„ = K„ - W„S„W„^ 
where v„ = y„ - y„ and W„ = 


Fig. 2. Non-augmented-state UKF O 













Initialise Set the prior mean //q and covariance Kq with a 
well-motivated estimate, along with parameters a, yS and 
y which determine the distribution of sigma-points a. 
Provide the estimated process and measurement 
covariances: and Zy. 


for n = l...n, - 1 do 

1) Augment the state mean and covariance K„_i 
with the the process noise and measurement noise 
means (zeros) and their respective covariances. 




0 
0 


, ^a,n = diag(K„_i,i:^,i:y). 


\ ^ / 

2) Calculate sigma points and weights 


N, + A _ 

X« = ^J(N, + A)Ka,n, 

= yu^ - y/(N. + d) K„,„, 


1 


2(N, + Ay 

for i - l...Nx, where A - (Nx + k), and Nx is the 
dimension of the augmented state vector. 

3) Transform sigma points using the forward map: 

f(x)+Xrf ^ 

x® = fa (x® where f^ (x^) = 


x^/ 

X,, 


, and 


X, Xd and Xy are the state, process-noise and 
measurement-noise parts of x^. 

4) Calculate the predicted mean and covariance: 

2N, 


fla,n = Yj 


?=0 


^a,n — (xi,i i^a,n) {^a^n Pfl.n) ■ 


i={) 

5) Apply the observation model to each transformed 
sigma point : 

y? = go (xl'l), where g„ (x^) = g (x) -H Xy. 

6) Calculate the predicted observation: 

IN-, 

i=0 

7) Calculate the innovation covariance: 


2N, 


s„ = 




/=0 


8) Calculate the cross covariance (for the 
non-augmented state): 

2N, 

i=0 

9) Perform a usual Kalman Filter update : 

f2n — f2n " 1 " 

K„ = K„ - W„S„W„^ 

where v„ = y„ - y„ and W„ = K^^S„'. 


end 


Fig. 3. Augmented-state UKF m, n 


Initialise Sample tip ‘particles’ from well-motivated 
prior: 

x®~P(x). 

Select initial weights VP® = 

Provide the estimated process and measurement 
covariances: Zd and Zy. 


for n - - 1 do 

1) Get next state sample: 

Propagate ‘particles’ through forward map, adding 
process noise sampled from relevant distribution 
(in this case, gaussian) 


Ji) 




2) Get log measurement probability: 

logp(y„|xW) = --vjz;^vi, where v,- = y„ - g(x®) 
and set weights 

= exp (log -H logp (y„|xW)). 


Normalise weights: VP„ 


(0 _ 


IP 


<0 


lljWli 


U) ■ 


3) Get effective number of particles: 

Weft = -^--y ■ 

if rieff < n,hr then 

Draw rip particles {.r,V^} from the current particle 


set 


|x®| with probabilities proportional to 


Reset weights - Ijiip. 


end 


end 


Fig. 4. Particle filter 


noise and measurement noise; the non-augmented state UKF 
samples the sigma-points twice, while the augmented-state 
UKF only once; also, in the non-augmented state UKF the 
process and measurement uncertainties are included in the 
calculation through the ‘augmented maps’ fa and ga (see Fig. 
1^ rather than added to the predicted covariance and innovation 
covariance as in the non-augmented UKF (see Fig. |^. 

The performance of these two algorithms is compared in 
Ref 0, which found that for the problems considered there 
the augmented-state UKF generally performs better. 

2) Particle filters: As well as the two varieties of UKF, we 
also test the performance of a particle filter. We use a sequen¬ 
tial importance resampling algorithm, first introduced in Ref. 
Cl. This algorithm has the advantage that as well as permitting 
a nonlinear forward map and observation function, it permits 
the prior distribution for the state vector to have any form, as 
well as the estimated noise distributions (although we do not 
need to exploit this flexibility for the current problem). The 
algorithm approximates a distribution of states by a sample of 
weighted states (‘particles’). It evolves each ‘particle’ through 
the forward map before performing a Bayesian update of the 
sample. Resampling prevents weights becoming concentrated 
in a small number of particles. 


















III. Results 


We ran the algorithms on measurements simulated using 
the following parameters/initial conditions (partly motivated 
by Ref. ISl): 6q - 0.2 rad, lOo - 0.2 rad/s, L - 3.5 m, mo = 1.7 
kg, mo = 5.5 kg/s, pA - 161.25 kg/m. We used - QA to 
generate process noise in the mass flow [via Eq. Q] at samples 
every 0.025s during 10s of filling time. Both UKF algorithms 
used a prior distribution for the state space with mean values 
equivalent to; 0o = 0.21 rad, wo = 0.15 rad/s, L - 2.5 m, 
niQ - 2.2 kg, mo = 5.36 kg/s, pA - 177.38 kg/m and variances; 
cr^ = 0.2 rad^ 


cr}, — 0.2 rad"^, erf ,,, 

ijA ’ logL/lm 

= 0.02, and cr^ 

logpA/lkg m 


= 0.2, cr} ,,, 

’ logm/lkg 


0.02, cr^ .,,, = 0.02, and cr^ ^= 0.2. This prior 

overlapped the ‘true’ state, but was not centred thereon. The 
particle filter sampled 1000 particles from a gaussian prior of 
the same mean and covariance as used by the UKF algorithms. 
For the estimated process noise covariance. Yd, we used for 
the log^m/lkg s ') variance component, and small but non¬ 
zero values for the other diagonal components. These non-zero 
values for the variables without expected process noise are 
necessary for the correct functioning of all the algorithms, and 
are usually justified in real applications by invoking the need to 
account for the difference between the (necessarily simplistic) 
model used and the laws governing the true system dynamics. 
In the simulated measurements, the noise covariance Yp - 
1.5 N^, and the estimated noise covariance supplied to the 
inference algorithms Yy - 2.5 

We ran all algorithms on the same set of simulated mea¬ 
surements illustrated in Fig. [TJb). The sequential estimates 
of mass are expressed as 90% prediction regions, and are 
shown in Fig. The most successful methods were the UKFs 
and particle-filter using the physics model [Figs, [^a), (c) and 
(e)]. These methods produced suitably narrow regions which 
contained the ‘true’ mass value. The UKFs which did not 
use a physics model provided wide 90% prediction intervals 
which contained the true mass value [Figs. |^b) and (d)]; the 
particle filter which did not use a physics model produced a 
narrow 90% prediction interval which often did not contain the 
true mass value [Fig. [^f)]. The particle filter without a model 
clearly attributed the oscillations in the measured vertical force 
component to oscillations in the mass flow, even though the 
magnitude of these oscillations was larger than would be 
expected from the stated process noise. We tested the non- 
augmented UKF for different realisations of the process and 
measurement noise, and it was found to produce consistently 
good estimates [see Fig. |^. 

As time f —» oo, the widths of the prediction intervals are 
limited by the effect of the finite process and measurement 
noise. By using idealised equipment, where the process and 
measurement noise is very small, it should be possible to 
obtain very accurate estimates of the mass. We demonstrate 
this by simulating a force dataset using negligible process 
noise in m. We ran all inference algorithms on this dataset 
using estimates of the process covariance cr^ with very small 
diagonal elements of order 10“* in each appropriate unit (and 
zero-valued off-diagonal elements). All inference methods 



Fig. 5. Sequential mass estimates (shaded regions) and true mass values 
[light (red) lines] as a function of time. Subplots (a) and (b) show results 
for the non-augmented-state UKF, (c) and (d) for the augmented-state UKF, 
and (e) and (f) for the particle filter. Results shown in (a), (c) and (e) are 
computed using the physics model and those in (b), (d) and (f) without. For 
the UKF results, grey regions represent the central 90% prediction regions 
for the mass, dotted lines represent the region boundaries and mean mass 
prediction; for the particle filter, grey regions represent masses between the 
5th and 95th percentile of the particle sample (indicated by dotted lines), the 
50th mass percentile is also given by a dotted line. 




Fig. 6. Sequential mass estimates (shaded regions) from for non-augmented- 
state UKF and true mass values [light (red) lines] as a function of time. Grey 
regions represent the central 90% prediction regions for the mass, dotted lines 
represent the region boundaries and mean mass prediction. Plots (a)-(d) show 
results obtained for simulated data with different realisations of the process 
and measurement noise. 


























Fig. 7. As Fig. but all with all inferences performed on simulated data 
with negligible process noise in m, using negligible estimated process noise. 

using the physics model produced estimates with very narrow 
prediction region [Fig. |^a), (c) and (e)]; however, the particle 
filter’s 90% prediction region did not contain the true mass 
value. For these simulations with no process noise [Figs.|7](b), 
(d) and (f)], the methods without a physics model produced 
poor estimates, either with wide prediction regions (in the case 
of the UKFs), or with narrow regions not always containing 
the true mass value (in the case of the particle filter). 

IV. Conclusions 

We have shown the benefit of including a physics model 
in sequential inference algorithms used in control-systems for 
automated weighing - in particular for weighing a suspended 
bag of powdered milk during the filling process. When used 
in conjunction with a UKF or particle filter algorithm, these 
methods are expected to produce narrow 90% prediction 
regions containing the true mass value, which would not be 
possible without the consideration of the dynamics generated 
by the physics model. Reference IQ found that the augmented- 
state UKF performed better than the non-augmented UKF 
in the problems considered there; however, for the problem 
considered in this paper, both UKF algorithms produce very 
similar output in similar running times. UKF algorithms are 
expected to be more useful than particle filters, since their 
running time is comparable to the filling times used in this 
paper. The particle filter’s running time is several times longer 
than these filling times, so in practice, particle filters would 
require slower filling rates and longer times between force 
measurements in order to inform a synchronous control sys¬ 
tem. Improved predictive performance of the particle filter in 


the cases where the prediction intervals fail to contain the true 
mass value would be expected if a larger number of particles 
were sampled. However, this would slow the running time 
even further. 

Consideration of systems with process noise in the density p 
and cross-sectional area A would provide straightforward and 
important extensions of this work. The density of milk would 
be expected to fluctuate as it settles during packing, and the 
cross-sectional area A would also be expected to vary with 
height within the bag. 
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